Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

The efect of the couple stress parameter and Prandtl number on the transient natural convection flow over a vertical cylinder

An analysis is performed to study transient free convective boundary layer flow of a couple stress fluid over a vertical cylinder, in the absence of body couples. The solution of the time-dependent non-linear and coupled governing equations is carried out with the aid of an unconditionally stable Crank–Nicolson type of numerical scheme. Numerical results for the steady-state velocity, temperature as well as the time histories of the skin-friction coefcient and Nusselt number are presented graphically and discussed. It is seen that for all flow variables as the couple stress control parameter, Co, is amplified, the time required for reaching the temporal maximum increases but the steady-state decreases.     

Issue Information

Three‐dimensional direct numerical simulation results of flow past a circular cylinder are influenced by numerical aspects, for example the spanwise domain length and the lateral boundary condition adopted for the simulation. It is found that inappropriate numerical set‐up, which restricts the development of intrinsic wake structure, leads to an over‐prediction of the onset point of the secondary wake instability (Re_cr). A best practice of the numerical set‐up is presented for the accurate prediction of Re_cr by direct numerical simulation while minimizing the computational cost. The cylinder span length should be chosen on the basis of the intrinsic wavelength of the wake structure to be simulated, whereas a long span length is not necessary. For the wake transitions above Re_cr, because the wake structures no longer follow particular wavelengths but become disordered and chaotic, a span length of more than 10 cylinder diameters (approximately three times the intrinsic wavelength) is recommended for the simulations to obtain wake structures and hydrodynamic forces that are not strongly restricted by the numerical set‐up. The performances of the periodic and symmetry lateral boundary conditions are compared and discussed. The symmetry boundary condition is recommended for predicting Re_cr, while the periodic boundary condition is recommended for simulating the wake structures above Re_cr. The general conclusions drawn through a circular cylinder are expected to be applicable to other bluff body configurations. Copyright © 2017 John Wiley & Sons, Ltd.     

Wave body interaction in 2D using smoothed particle hydrodynamics (SPH) with variable particle mass

Wave interaction with bodies is an important practical application for smoothed particle hydrodynamics (SPH) which in principle applies to steep and breaking waves without special treatment. However, few detailed tests have been undertaken even with small amplitude waves. In order to reduce computer time a variable particle mass distribution is tested here with fine resolution near the body and coarse resolution further away, while maintaining a uniform kernel size. We consider two well‐defined test cases, in two dimensions, of waves generated by a heaving semi‐immersed cylinder and progressive waves interacting with a fixed cylinder. But first, still water with hydrostatic pressure is tested. The open‐source code SPHysics ( http://www.sphysics.org )^{§Update made here after initial online publication.} is used with a Riemann solver in an Arbitrary Lagrangian–Eulerian formulation. For the heaving cylinder, SPH results for far field wave amplitude and cylinder force show good agreement with the data of Yu and Ursell (J. Fluid Mech. 1961; 11 :529–551). For wave loading on a half‐submerged cylinder the agreement with the experimental data of Dixon et al. (J. Waterway Port Coastal Ocean Div. 1979; 105 :421–438) for the root mean square force is within 2%. For more submerged cases, the results show some discrepancy, but this was also found with other modelling approaches. The sensitivity of results to the value of the slope limiter used in the MUSCL‐based Riemann solver is demonstrated. The variable mass distribution leads to a computer run speedup of nearly 200% in these cases. Copyright © 2011 John Wiley & Sons, Ltd.     

A ghost-cell immersed boundary method for large-eddy simulations of compressible turbulent flows

A methodology to perform a ghost-cell-based immersed boundary method (GCIBM) is presented for simulating compressible turbulent flows around complex geometries. In this method, the boundary condition on the immersed boundary is enforced through the use of ‘ghost cells’ that are located inside the solid body. The computations of variables on these ghost cells are achieved using linear interpolation schemes. The validity and applicability of the proposed method is verified using a three-dimensional (3D) flow over a circular cylinder, and a large-eddy simulation of fully developed 3D turbulent flow in a channel with a wavy surface. The results agree well with the previous numerical and experimental results, given that the grid resolution is reasonably fine. To demonstrate the capability of the method for higher Mach numbers, supersonic turbulent flow over a circular cylinder is presented. While more work still needs to be done to demonstrate higher robustness and accuracy, the present work provides interesting insights using the GCIBM for the compressible flows.     

A ghost fluid Lattice Boltzmann method for complex geometries

A ghost fluid Lattice Boltzmann method (GF‐LBM) is developed in this study to represent complex boundaries in Lattice Boltzmann simulations of fluid flows. Velocity and density values at the ghost points are extrapolated from the fluid interior and domain boundary via obtaining image points along the boundary normal inside the fluid domain. A general bilinear interpolation algorithm is used to obtain values at image points which are then extrapolated to ghost nodes thus satisfying hydrodynamic boundary conditions. The method ensures no‐penetration and no‐slip conditions at the boundaries. Equilibrium distribution functions at the ghost points are computed using the extrapolated values of the hydrodynamic variables, while non‐equilibrium distribution functions are extrapolated from the interior nodes. The method developed is general, and is capable of prescribing Dirichlet as well as Neumann boundary conditions for pressure and velocity. Consistency and second‐order accuracy of the method are established by running three test problems including cylindrical Couette flow, flow between eccentric rotating cylinders and flow over a cylinder in a confined channel. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of adjustable permanent magnetic field sources

A permanent magnet assembly in which the flux density can be altered by a mechanical operation is often significantly smaller than comparable electromagnets and also requires no electrical power to operate. In this paper five permanent magnet designs in which the magnetic flux density can be altered are analyzed using numerical simulations, and compared based on the generated magnetic flux density in a sample volume and the amount of magnet material used. The designs are the concentric Halbach cylinder, the two half Halbach cylinders, the two linear Halbach arrays and the four and six rod mangle. The concentric Halbach cylinder design is found to be the best performing design, i.e. the design that provides the most magnetic flux density using the least amount of magnet material. A concentric Halbach cylinder has been constructed and the magnetic flux density, the homogeneity and the direction of the magnetic field are measured and compared with numerical simulation and a good agrement is found.     

Limit cycles and invariant cylinders for a class of continuous and discontinuous vector field in dimention 2n

The subject of this paper concerns with the bifurcation of limit cycles and invariant cylinders from a global center of a linear differential system in dimension 2n perturbed inside a class of continuous and discontinuous piecewise linear differential systems. Our main results show that at most one limit cycle and at most one invariant cylinder can bifurcate using the expansion of the displacement function up to first order with respect to a small parameter. This upper bound is reached. For proving these results we use the averaging theory in a form where the differentiability of the system is not needed.     

Influence of power-law index on transitional Reynolds numbers for flow over a semi-circular cylinder

In this work, the governing partial differential equations (continuity and Cauchy’s momentum equations) describing the flow of power-law type non-Newtonian fluids across a semi-circular cylinder (oriented with its curved surface in the upstream direction) have been solved numerically. In particular, consideration has been given to the delineation of the critical Reynolds numbers denoting the onset of flow separation from the surface of the cylinder and the onset of the laminar vortex shedding regime. This information is germane to establish the scaling of the macroscopic characteristics like drag coefficient and Strouhal number on the governing parameters, namely, Reynolds number and power-law index. The present results clearly suggest that the transitional Reynolds numbers show a strong dependence on the type (shear-thinning and shear-thickening) of fluid behavior as well as on the severity of the shear-dependence of the viscosity. With reference to the behavior seen in Newtonian fluids, the flow remains not only attached to the surface up to higher Reynolds numbers, but shear-thinning behavior also delays the onset of the laminar vortex shedding regime. As expected, shear-thickening fluids, of course, display the opposite characteristics.     

Improvements on open and traction boundary conditions for Navier–Stokes time-splitting methods

We present in this paper a numerical scheme for incompressible Navier–Stokes equations with open and traction boundary conditions, in the framework of pressure-correction methods. A new way to enforce this type of boundary condition is proposed and provides higher pressure and velocity convergence rates in space and time than found in the present state of the art. We illustrate this result by computing some numerical and physical tests. In particular, we establish reference solutions of a laminar flow in a geometry where a bifurcation takes place and of the unsteady flow around a square cylinder.